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Abstract 

We reexamine non-Einsteinian effects observable in the orbital motion of 
low-orbit artificial Earth satellites. The motivations for doing so are twofold: 
(i) recent theoretical studies suggest that the correct theory of gravity might 
contain a scalar contribution which has been reduced to a small value by the 
effect of the cosmological expansion; (ii) presently developed space technolo- 
gies should soon give access to a new generation of satellites endowed with 
drag-free systems and tracked in three dimensions at the centimeter level. Our 
analysis suggests that such data could measure two independent combinations 
of the Eddington parameters [i = (5 — 1 and 7 = 7 — ! at the 10~^ level and 
probe the time variability of Newton's "constant" at the G/G ~ 10~^^yr~^ 
level. These tests would provide well-needed complements to the results of the 
Lunar Laser Ranging experiment, and of the presently planned experiments 
aiming at measuring 7. In view of the strong demands they make on the 
level of non-gravitational perturbations, these tests might require a dedicated 
mission consisting of an optimized passive drag-free satellite. 

PACS number(s): 04.80.Cc, 04.25.Nx, 95.40.+S 
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I. INTRODUCTION 



Though relativistic gravity has already been probed by quite a few experimental tests, 
with the conclusion that General Relativity is in agreement with all existing experiments 
(see for reviews), recent theoretical developments in tensor-scalar cosmological models 
suggest that the correct theory of gravity might differ from Einstein's theory by containing 
a scalar contribution whose present magnitude has been naturally driven to a small value 
by the cosmological expansion |^,^. These theoretical results provide a new motivation for 
trying to improve the existing tests of tensor-scalar theories of gravity. 

In this paper we reexamine the use of artificial (Earth) satellites to perform orbital tests 
of non-Einsteinian gravity. We shall consider here only the most conservative deviations from 
Einstein theory, i.e. boost-invariant effects associated to the exchange of scalar excitations. 
[We refer to our recent work [Q for a discussion of the use of satellite data as probes of 
possible preferred-frame effects in relativistic gravity.] Although several of the Einsteinian 
or non-Einsteinian effects in the motion of artificial satellites have been abundantly discussed 
in the literature , it seems that there exists no treatment of the problem which is both 

systematic and observation-oriented. There exists a systematic analysis of non-Einsteinian 
effects in the motion of the Moon [^j, but we shall see that the two problems are quite 
dissimilar both because of differences in physical parameters (aMoon/-REarth ~ 60 compared 
to dsateiiiteZ-REarth ~ a fcw) and of differences in observational techniques (radial ranging for 
the Moon versus three-dimensional tracking for satellites). The aim of the present paper is 
to provide a systematic analysis of non-Einsteinian effects possibly observable in the orbital 
data of drag-free satellites whose three-dimensional motion is tracked (either by Laser or by 
Global Positioning System (GPS) techniques) at the centimeter level. The Gravity Probe B 
(GPB) satellite should soon provide the first example of such satellites. [The existing Laser 
Geodynamical Satellites (LAGEOS) are tracked with centimeter precision, but the small 
level of residual drag acting on their motion is probably too large for them to be useful 
probes of the effects we discuss below.] 

Note that, contrary to some proposals, e.g. [0, we shall not consider explicitly here the 
purely Einsteinian relativistic effects. We consider that the ensemble of existing positive 
tests of General Relativity [||,|| has already probed most of the qualitatively crucial fea- 
tures of Einstein's theory. We shall therefore assume that the Einsteinian post-Newtonian 
contributions to the satellite equations of motion are separately included in their entirety 
(see ^^), and focus on the quantitatively most promising effects for measuring possible 
non-Einsteinian weak-field deviations. 

In the post-Newtonian limit, generic tensor-scalar theories of gravity predict that the 
motion of spherical, non rotating bodies is given by the Lagrangian [III 



L-,,,rs = -T.mAc\l-.l/c'y/' 
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-\ E (1 + 2/3)^!?^!^, (1.1) 

where tab = \^a — ^b\, ^ab = {^a — ^B)/rAB, ^a = d:KA/dt, and where the summation 
symbols denote muhiple sums over the various independent body labels A, 5,. .. with the 
only exclusion of the would be infinite contributions. 

In Eq. ( |1.1| ) the effective coupling constant Gab for the gravitational interaction of bodies 
A and B is of the formQ 

Gab = G[1 + 6a + 6b] (1.2) 
where Sa is generically the sum of two physically independent contributions 

_ pgrav 

Sa = Sa + V^, (1-3) 
rriAC^ 

with 

r] = Ap-^. (1.4) 

The first contribution 6a comes from the fact that the physically best motivated tensor- 
scalar theories violate the weak equivalence principle 0] , while the contribution proportional 
to the gravitational self-energy 

Ef''^ = -{G/2) [ [ d'xd'x' p{i^)p{id)l\^-^\ 

J A J A 

derives from the fact that all the theories containing a variable local gravitational "constant" 
violate the strong equivalence principle |1^ . 



The dimensionless numbers /?, 7 parametrize weak-field deviations from General Rel- 
ativity and are proportional to the strength (relative to the usual tensor contribution) of 



the coupling to matter of the scalar contribution to gravity [jTT|. They are related to the 
usual post-Newtonian (Eddington) parameters simply by/5 = /3 — 1,7 = 7 — l(so that 
Tj = 4/3 — 7 — 3). Besides the effects parametrized by /5, 7 and the 5's, another deviation from 
General Relativity generically predicted by tensor-scalar theories, and implicitly contained 
in the Lagrangian ( |1 . 1| ) , is a possible time variability of the gravitational coupling strength 
G, as well as (in generic tensor-scalar theories 0]) a time variability of the massesQmA- 



^In Eq. (|1.2| ), the basic gravitational constant G is the sum of the contribution due to the tensor 
interaction, and of a mean contribution due to the scalar interaction. One assumes that the small 

— —2 

composition-dependent effects 6a ^ ^, and neglects terms of order 6 . 

^The other parameters (3, 7, 6a are generically also slowly changing. However, this leads to 
higher-order, unobservably small effects. 
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The best present (one sigma) observational limits on the parameters of the Lagrangian 
( |1.1|) measuring deviations from Einstein's theory are: I7I < 2 x 10"'^ |jl3|, |/?| < 3 x 10~^ 
[assuming a Sun's J2 ~ 2 x 10"''') ||14|| , ^Moon ~ ^Earth = (-2.7 ± 6.2) X 10-13 ig, and 



\G/G\ ^ 10 iiyr ^ [|T4| (we are using the fact that combined geophysical and astronomical 
data set a bound of order lO^^^yr^^ on the possible time variability of the masses |16[)- Note 



that if one were to assume the exact validity of the weak equivalence principle, i.e. 5a = 0, 
the Lunar Laser Ranging result on ^Moon — ^Earth would imply r] = (—0.6 ± 1.4) x 10-^, from 
which — using the Viking limit — one would deduce (3 = (—1 ± 5) x 10"^ [jl5[. However, 
in the present work we have in mind the general class of tensor-scalar theories which do not 
respect the weak equivalence principle (such as in Ref. [Q). As the present observational 
limits on the 5^'s are not better than 3 x 10^^^ |l^, we must consider the S^s in Eq. ( [L . 1|) 
as being independent from (3 and 7. 

The object of the present paper is to point out that the centimeter-level tracking of 
drag-free artificial Earth satellites is a promising tool for measuring both (3 and 7 at the 
10—^ level, and G/G at the lO^^^yr^^ level. We shall show that there are two independent 
relativistic effects which build up over many orbital periods and lead to displacements ~ 
0(/3,7) X lO^cm for integration times of the order of 6 months. One of these two relativistic 
effects is non null and is the well known perigee advance associated with the (Eddington- 
modified) Schwarzschild field of the Earth which has been abundantly discussed in the 
literature |l|J^,|,|l^. The other one is a null relativistic effect whose origin is known [|l],3, 
but whose importance, in the context of low orbit satellites, for testing the combination 
7] = 4(3 — ^ has (as far as we are aware) not been clearly realized before. [See further 
discussion below. 1 



II. LAGRANGIAN AND EQUATIONS FOR MOTION FOR SATELLITE 

MOTION. 

The Lagrangian (|1 . 1| ) describes the global dynamics of the bodies of the solar system 
( "barycentric" reference frame). We wish to pass to a "geocentric" reference frame. A 
rigorous and elegant way of effecting such a transformation in the framework of General 



Relativity has been recently discussed in Ref. [T^. However, we are here working within a 
non-general-relativistic framework. The safest way of deriving the deviations from a general 
relativistic description of the geocentric motion of a satellite is to start from the /? — 7 — 5- 
modified barycentric equations of motion of the Earth and the satellite and to take their 
difference. We restrict our attention to the effects associated with the satellite (label 1), the 
Earth (label 2) and the Sun (label 3), and introduce the notation^ r = r n = X12 = Xi — X2, 
V = r, DN = X32 = X3 - X2, V = X32 = V3 - V2, m = m2 and M = m^. 

The general relativistic equations of motion exhibit some remarkable features called "ef- 



facement properties" in Ref. |19]. There is an effacement of the internal structures of the 



gravitating bodies in the sense that, to a high accuracy, one can write the equations of 



We do not need a special notation for the mass of the satellite as it disappears from the final 
equations. 
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motion in terms of only some "centers of mass" and some "masses" . There is also an efface- 
ment of the external universe in the dynamics of a local gravitating system. For instance, 
after introducing some special geocentric coordinate frame, the influence of the Sun on the 
dynamics of an Earth satellite is effaced up to terms which decrease like when D oo. 
Both types of effacement properties are violated in tensor-scalar theories. The violation of 
the first type of effacement properties shows up clearly in the Lagrangian ( [L.l|) through the 
composition dependence of the effective coupling constant Gab- The violation of the second 
type of effacement properties shows up in the fact that the geocentric equations of motion of 
a satellite will contain terms proportional to GM/c^D and GM/c^D"^ (with coefficients lin- 
ear in (3, 7 and 6) which cannot be eliminated by a coordinate transformation. In this paper, 
we shall be especially interested in the leading (oc GM/c^D) and subleading (oc GM/c^D"^) 
violations of the effacement of external influences appearing in non-Einsteinian theories. We 
shall neglect the terms containing a factor GM/c^D^ which are, roughly speaking, non- 
Einsteinian relativistic corrections to the Newtonian tidal forces and are therefore extremely 
small0. 

It is convenient to simplify the geocentric equations of motion by rescaling the local 
spatial coordinates according toQ 



1 + 7 



GM \ 



,old 



(2.1) 



After this rescaling, the geometric equations of motion have the form 



Agr + B 



(2.2) 



with 



^GR 



Gm GM ^ ^, 

n - ^TT^ [r - 3(r ■ N)N] 



(2.3) 



denoting the general relativistic acceleration term (see [10| for the full expression of Aqr 
including a relativistic treatment of tidal forces), and 



B = Bo + Bi + B2 + 



GM f(3_ 7 -^ 

JJZ L2 ^ C2 ^ I 



(2.4) 



denoting the expansion of the non-Einsteinian acceleration terms in successive powers of 



^Compared to the leading Newtonian acceleration, they are typically smaller by a factor ^ 
^{V/cf Mr^/mD^ < 5 x IQ-^^ {a / R^^rthf . 

^This rescaling is the 7-dependent part of a standard rescaling in the PPN literature [|l|,§. Its 
effect is to efface the term proportional to ^GM/c^D in the spatial metric gij. 
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Br 



Gm 



n 



Gm V 
2(/5 + 7)-^-7- 



(5i - (52 



27 Gm , 
+ — — (n ■ v)v 



Bi 
B2 



Gm Gm 
(4/j-7) „ n 



GM 



N 



"-2^^ + 5,-^2 



(2.5) 
(2.6) 



c^r 



^,-5 .GMGm. 
+2(/3 + 7)-— ^(n-N)n 



+2 



^GM 
^ GM 



^2 (^2j. 

v^N - 2(v ■ N)v 
(N X V) X V . 



(2.7) 



To illustrate the non-Einsteinian terms proportional to GM/ we can quote its numer- 
ically dominant contribution: 



B3 = -(5i + 53)^[r - 3(r ■ N)N] + 



(2.8) 



Let us also note that we shall neglect in this work the non-Einsteinian acceleration terms 
depending on the higher (mass and spin) multipole moments of the Earth. [We assume 
however that the Einsteinian effects of these is fully taken into account in Aqr.] For instance, 
the spin-orbit interaction adds a term of the form B5 

(GS2 X r/r^) where S2 is the spin angular momentum of the Earth, so that 



V X H- to B. Here Hy = — x 



B< 



7 



G 

3 



l^v X S2 — 3(n • S2)v X nj 



(2.9) 



The spin-orbit acceleration B5 is smaller than the velocity-dependent terms in Bq, Eq. ( p. 51) , 
by a factor S2/{mrv) ~ 0.02 (a/i?Earth)~^''^- The other multipole contributions will contain 
even smaller factors. They can be therefore safely neglected with respect to the leading 7- 
and /3-dependent terms kept in Eqs. (|2.5|)-(|2r7|). 

The last term on the right-hand side of Eq. ( p.7|) has the form of a Coriolis force, 20-x v, 
with 



GM 

"^ = ^^^^^- 

As first discovered by De Sitter a similar relativistic Coriolis term, 2f2GR x v, with 

3 GM 



GR 



N X V 



(2.10) 



(2.11) 



is present in Einstein's theory if one works in a geocentric frame which has a fixed orientation 
with respect to the barycentric frame. Alternatively, one can efface away the total relativistic 
Coriolis force from the local equations of motion by working in a geocentric frame which 
rotates with angular velocity 
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tot 



with respect to the barycentric frame 
the one given here.] Let us note that the rotation fitot 



[See 



2_\ 

1 + ^7 ) f^GR 



(2.12) 



T0|] for a more accurate treatment of flcR than 
"precession of the local Earth inertial 



frame") is acting upon all the satellite orbits (including the Moon), as well as upon the solid 
Earth. Therefore Otot essentially disappears in the position measurements effected by means 
of GPS or Laser techniques. It appears however when connecting local measurements to 
global ones, e.g. through the use of Very Large Baseline Interferometry (VLBI) techniques. 
It appears also in the Lunar Laser Ranging (LLR) data because the internal dynamics of the 
Earth-Moon system is strongly influenced by the Sun. The measurement of fitot through 
LLR data is now reaching the 1% level (more precisely (2/3)7 = —0.3 ± 0.9% [^). As 
the effects of fitot on low satellite orbits are more difficult to observe^ this does not look 
as a promising way of getting improved tests of 7. Therefore, in the following, we shall 
concentrate on the other non-Einsteinian effects B in Eq. ( ^.21 ). In other words, we shall 
study the satellite motion in a geocentric frame precessing with the velocity ( |2.12|) , so that 
the last term in Eq. ( |2.7|) is absent. 

In such a dynamically non rotating geocentric frame one can verify that the equations 
of motion of a satellite derive from a Lagrangian of the form 



+ R 



(2.13) 



1 9 Gm 
-v^ H 

2 r 



GM 
2D^ 



[r^ - 3(N ■ r) 



+ other Newtonian and post-Newtonian contributions. 



(2.14) 



R = R0 + R1 + R2 + O 



GM 



(2.15) 



Rq 

Ri 
R2 



— - V — Gm 

51 + ^2 + 7^-/3 



c^r 



^_ GM Gm 

c U r 
^ Gm - 



7^-2/3 



+ 5i 



Gm 



-\GM.^^ , 
52N^(N-r) 



(2.16) 
(2.17) 
(2.18) 



^They correspond to spatial displacements Ototo = (1 + |7)(a/-REarth) x 59 cm yr 



7 



III. NON-EINSTEINIAN EFFECTS IN THE MOTION OF SATELLITES 



As we are discussing small deviations from a general relativistic motion, we can discuss 
independently the perturbations depending upon different powers of in Eq. (|2.15|) and 
superpose linearly their effects. 

The terms which are independent from the presence of the Sun, Eqs. ( p.5|) and ( p.l6|) , are 
well known. Apart from a constant renormalization of G ^ ^effective = G[l + 5i + ^2] (whose 
possible dependence upon the satellite's composition is too small to be of observational sig- 
nificance), they correspond to a geodesic motion in an Eddington-modified Schwarzschild 
metric. The corresponding solution can be written down explicitly (see e.g. [^] and Ap- 
pendix A. 2. 3 of for a simple form of the explicit solution at the post-Newtonian accuracy). 
Short-period perturbations correspond to spatial displacements of order {j3 + ^)Gm/c'^, i.e. 
smaller than about 10~'^cm given the existing limits on (3 or 7. This is too small to be of 
observational significance. Finally, the only observationally significant effect coming from 
Bq is the secular advance of the perigee 

where a denotes the semi-major axis of the orbit, e its eccentricity and n = 2tt/P ~ 
{Gm/ 0?)^/"^ the orbital frequency [The index zero on the left-hand side of (3.1) refers to 
the cause Bq of this effect]. In terms of the ratio a = a/R where R = 6.371 x 10®cm is the 
Earth radius, one findsQ 

- a~3/2 / t \ 

aSouj ~ (27 - f3) X X 1.74 X lO^cm . (3.2) 

1 — \ 1 yr y 



The terms which are proportional to GM/D, Eqs. ( p.6|) and ( p. 171) are, in principle, well 



known. They can be simply interpreted as associated with a renormalization of the local 
(geocentric) value of G due to the proximity of the Sun: G — > Gioc = (1 — rjGM / c^D)G 
The existence of this (time-dependent) renormalization was taken into account in the 
treatment of non-Einsteinian effects in Lunar Laser Ranging but was found to lead to a 
very small perturbation of the (directly observed) radial distance (about 140 times smaller 
than the effects associated with the violation of the equivalence principle, — ^2 7^ 
in Eq. ( p. 7] )). This conclusion is correct in the case of the Moon, but the main point of 
the present paper is to emphasize that the situation is very different in the case of low 
Earth satellites for two different reasons: (a) their proximity to the Earth enhances the 
perturbation, and (b) the fact that their motion can be tracked in three dimensions gives 
access to the longitudinal displacement which happens to be enhanced by a large factor 
oc n/n2, where n2 = {GM / a^y/"^ denotes the orbital frequency of the Earth around the Sun. 



^Numerically P = a^/^ x 1.406 h, Gm/c^ = 0.4435 cm, Gm/Rc^ = 6.96 x 10" 



8 



It is possible to solve for the effects of Ri, Eq. (|2.17| ), by the following method. Let us 
introduce new space and time variable^ 

GM\ 

l-r^^jr, (3.3a) 

t' = t- e'^ = t-2r] [ dt ^ . (3.3b) 

The difference between the Lagrangian L' describing the motion in the new variables and 
the original Lagrangian L is given by ||2l| 



L'-L=^-^{e-e\) + ^Q{e) , (3.4) 
or dt 

where 5L/5r = —{r + Gmn/r^) + 0{GM/D'^) + 0{l/c^) and where dQ/dt is some total 
time derivative. Using the identities 

^ n \ , Gm 1 c/^r^ 
r + Gm — ] = -v^ + + 



r^/ r 2 dt"^ 



^ n\ d (\ Gm\ 



r^/ dt \2 r 

and operating by parts one finds 



+ 0{^e]+0[-) . (3.6) 



GM \ ^( e 
6 

The first term on the right-hand side of (|3.6| ) cancels the contribution i?i to L, while the sec- 
ond term is easily seen to be of order rjiV / c^GMr'^ / D'^ , i.e. of order of the non-Einsteinian 
relativistic corrections to tidal effects (see B3, Eq. (|2.8|)) that we are neglecting in this work. 

The conclusion is that the effect of Ri is equivalent to performing the transformation 
( ^.31 ) on the solution r'(t') of the equations of motion in which Ri (and Bi) have been set 
equal to zero. Note that this method can be applied either to the full term GM/c^D or only 
to its time-varying piece due to the eccentricity 62 of the Earth orbit around the Sun (we 
neglect (62)^ for simplicity): 

GM GM GM GM , , , , 

where ^2 is the time of passage of the Earth at its perihelion. The time-independent renor- 
malization of G associated with GM/c^a2 is locally unobservable. The time-dependent 



^We do not consider r' and t' as new coordinates in space-time but only as auxiliary variables 
which are helpful to solve the equations of motion. Note that t' would have bad properties as a 
time coordinate as it would introduce a term linear in r]GM/c^ D in the metric component ^qo- 



9 



perturbations in the spatial positions, Eq. ( p.3a|) , are of order O.lrja cm, i.e. too small to be 
observable. Finally, the only observationally significant effect coming from Bi is the change, 
given by Eq. ( p.3b| ), in the time according to which the satellite runs on its orbit. If we 
introduce the mean anomaly i, i.e. an angle connected with the position on the orbit which, 
in the unperturbed motion, varies linearly in time, we shall have in the perturbed motion 



nt + const = n [t — 2r] j at j + const , (3.8) 



OT i = nt + Si£ with 



GM n . . , , 

Oil = —27] 62 -7: sin n2(t — 12) + const . (3.9) 

c^a2 n2 

Using 62 = 1.673 x 10"^ and GM/c^a2 = 9.87 x 10"^, the perturbation corresponds to 
an along-track displacement 

aSii = -(J3 -l^)a-^^'^sm 1 27r^^ | x 5.25 x 10^ cm . (3.10) 
4 \ 1 yr / 

The yearly modulation of Gioc = (1 ^ rjGM /c^D)G has other consequences which, in 
principle, enter satellite data through the motion of, e.g.. Laser stations on the ground 
used to track a satellite. First, the radius of the Earth will undergo a corresponding yearly 
modulation. Using dhiR/ d\n.G ~ —0.1 ||22|, the amplitude of this modulation is 

5li?^-0.1i?^^0.1i?62 ^ b - 77) COS 712 {t-t2) 

Gioc c^a2 V 4 / 

~^ - ^7^ cos (^Tx ^ ^-21 ^ cm , (3.11) 

which is negligibly small. Second, this yearly breathing of the Earth will entail a cor- 
responding modulation of its angular velocity, say VL2. Using 80.2/^2 = —8X2IT2 and 
5X2/T2 — —0.VJ5G/G 1^ where I2 denotes the Earth inertia moment, this leads, by 
integration, to a yearly modulation of the angle of rotation of the Earth and therefore a 
corresponding longitudinal displacement of ground stations: 

^2 4GM 
— ^2 — — . . 

712 6^02 V 4 

^7) sin (^Ti ^ ^^-^ ■ '^^■^^^ 

Again, this effect is negligible compared to the orbital effect (|3.10| ). [It would moreover 
probably be impossible to distinguish from other yearly modulations in the rotation of the 
Earth.] 

It remains to study the orbital effects of the terms proportional to GM/c^D"^ in the 
equations of motion. One checks that the short-period effects associated with B2, Eq. (p. 71), 
are too small to be observationally relevant. One can get the secular effects of the orbital 



R5i^^-0.l7R — e2— — -7)sinn2(t-t2) 



10 



elements caused by B2 by averaging over one orbital period the time derivatives of the 
energy, the angular momentum and the Lagrange-Laplace (-Runge-Lenz) vector. One finds 
that {da/dt)2 = and 

{de/dt)2 = f2 X € , (3.13a) 
{di/dt)2 = f2 X e , (3.13b) 

where e denotes the eccentricity vector of the orbit (i.e. a vector of magnitude e directed 
toward the perigee), £ = ^/l — c, with c a unit vector along the orbital angular momentum, 
and where the average "forcing term" caused by B2 is 



GM 

3 ^ \ na 



^2 = k—^N, (3.14a) 



..^(.,-..J--(4/3-T). (3.14b) 



In the general case of eccentric orbits, Eq. (^.13b|) shows that B2 will produce variations of 



the inclination I and the node Q to the orbit of the satellite. Here, we shall concentrate 
on the more usual case where the eccentricity is small. Then, Eq. ( |3.13a| ) with £ ~ c gives 



the contribution of B2 to the variation of e. We must add the Newtonian effects of the 
quadrupole moment of the Earth ( J2 = 1.08263 x 10~^). Let us introduce the basis (a, b) of 
the orbital plane where a is directed toward the ascending node and a x b = c. Because of 
the Newtonian effects of J2 this basis is rotating around the Earth's polar axis with angular 
velocity Cl ~ —^n J2d~^ cos J (see e.g. 0). Let d' /dt denote a time derivative in the rotating 
frame (a, b, c). We get the following equation for the secular time evolution of the two 
independent components of e with respect to the vectors (a, b) 

d'e/dt = c X [unG- i^{t)] , (3.15) 

where = (f2 ■ a)a+ (f2 ■ b)b denotes the projection of £2 onto the orbital plane, and where 

3 nJ2 , . r ■ 2 T\ 

^N = (4-5 sm /) 

4 

~ (4 - 5 sin^ /)a~^/2 x 27r/(0.20 yr) (3.16) 



is the Newtonian perigee advance due to the Earth's quadrupole moment Q]. Eq. (|3.15 



constitutes an inhomogeneous linear differential equation in e = ecosa;a + e sin cub. The 
time dependence of the forcing term f^(t) comes from a combination of the yearly variation 
of f2, Eq. ( p.l4a|) , with the rotation of the orbital plane with angular velocity 17 around the 



Earth's polar axis. Following our recent work 0|, we can easily solve for e(i(:) by decomposing 
f_L in a sum of constant-norm vectors fa rotating with constant angular frequencies cUq in the 
(a, b)-plane: f_L = Eafo,. Alternatively, one can work with complex numbers, xa + ?/b — >■ 
z = X \ iy, and decompose f_L ■ (a + ib) into harmonic components oc exp(ia;Q,t). Then the 
general solution of ( |3.15| ) has the form 

e(t) = e;v(t) + Eea(i) , (3.17) 
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where 



(3.18) 



and where e^(t) is a constant-norm vector rotating in the (a, b)-plane with angular velocity 
un (homogeneous solution of (3.15)). In terms of complex numbers, z = e ■ {a + ih) = 
e expiuj, this general solution reads 



[e(t) exp iuj{t)]2 = Cn ^W^i^Nt + ctn) 



+k 



GM 
4 



^1 + COS/2)(l + cos/) 



cutv + ^ — n2 
1 + cos 12) (1 — cos/) 



exp +i{n2t — Q) 



cjjv — f2 + n2 

'1 — COS/2)(l + cos/) 



+ 



ujn + ^ + n2 

[1 — COS/2)(l — cos/) 



+ 



ujn' — ^ — n2 
2 sin I2 sin / 



exp —i{n2t — Q) 
exp —i{n2t + VL) 
exp +i{n2t + VL) 



— n2 
2 sin I2 sin / 



UJn + "^2 



exp +1 n2t 



exp —i n2t 



(3.19) 



where I2 = 23.45° is the inclination of the apparent orbit of the Sun with respect to the 
Earth equatorial plane, and where we have neglected the eccentricity 62 of the Earth's orbit. 

We see clearly from Eq. ( |3.19| ) that the effects of B2 can be enhanced by small divisors 
if the semi-major axis and the inclination of the satellite orbit (on which both un and 
Cl depend) are appropriately chosen. The locus of these resonances in the (a, /) plane is 
the same as for the perturbations due to the orbital velocity of the Earth around the Sun 
in non-boost-invariant gravity theories. It will be found in Fig. 3 of Note that these 
resonances exist only if a is smaller than 2.65/?. Because of the smallness of I2, the dominant 
contributions in the square brackets on the right-hand side of Eq. ( p.l9| ) are the first two ones 
corresponding to the frequencies ±(n2 — Q) (see the solid lines in Fig. 3 of 0). However, the 
price to pay to take advantage of these resonances is to dispose of a long time of observation. 
In this work we have in mind using tracking data from missions dedicated to other purposes 
(e.g. GPB or STEP) which may last only for a year or so. Then, the best situation is when 
uun ± (f^ — ■"-2) is of order 27r/l yr. In such a case, we see from Eq. (p.l9|) that B2 will cause 
spatial displacements with periods ~ 1 yr and with amplitudes 



0^26 ~ ak 



GM 1 ± cos / 
2/^2 Co^±{Q- n2) 



[1.80 X 10^ a'/\6i - 62) - 1.67 a^/\f3 - ^7)] 
27r/l yr 

X ( 1 ± cos /) : cm , 

CUAT ± (fi - 77-2) 



(3.20) 
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where 51-62^61-62 + 1.86 x 10-9(/3 - 7/4), cf. Eq. For low Earth sateUites (a ~ 1) 

this is much smaller than the effects of Bq and Bi discussed above, and will therefore not 
contribute much to the measurability of /?, 7 and 6. By comparing however the negative 
power law dependences upon a of Eqs. ( |3.2| ) and ( |3.10[ ) with the corresponding positive power 



dependences in Eq. (|3.2CI| ) we see in retrospect why the most sensitive non-Einsteinian test 



in the motion of the Moon (d ~ 60) can be the violation of the equivalence principle (term 
<x6i-62 mk^g). 

Finally, let us consider for completeness the orbital effects of a slow time variation of 
G and/or the masses. Going back to the level of the original Lagrangian ( |1.1| ) we must 
consider that the mass of the satellite rris, the mass of the Earth m and G can vary in the 
Newtonian-approximation Lagrangian for the geocentric motion: 




G{t)m,{t)m{t) _ 



We can solve for the dynamics of Lat by a simple generalization of the method used 
above for the yearly variation of Gioc. 

Indeed, let us consider more generally a perturbed Lagrangian of the type L = Lq + Li 
with Lo = Iv^ + Gm/r and Li = ^a{t)Y'^ + b{t)Gm/r where Gm is constant and where a{t), 
b{t) ^ 1. Let us introduce the following new space and time variables 

r' = (l + a + 6)r, (3.22a) 
t' = t + J dt{a + 26) . (3.22b) 

By using Eqs. (|3.4| ) and (|0| ) we find that 

We can apply this result to the Lagrangian (|3.21| ) decomposed as Ln = rUslLo + Li\ by 
writing ms{t) = 771^(1 + trhs/ms)-, m{t) = m(l + tm/m), G{t) = G{1 + tG/G) (where m^, m 
and G denote some constants). [Here, one is making the assumption that the characteristic 
time scale T of variation of G{t), ms{t), m{t) is very large so that one can work to first 
order in T^^.] This yields a = trhs/ms and h = t{G/G + m/m + Ths/mg), and therefore 
{d? /dt^){a + 6) = in Eq. (|3.23|) . Finally, we get that the motion in terms of the variables 



( G fn ^rhs 
1+t - + — + 2— 
V G m m. 



(3.24a) 




2^ + 2^ + 3^] , (3.24b) 
G m msj 

is described by the unperturbed Lagrangian Lq. 

Knowing already that \G/G\ < IQ^^^ yi~^ and \rri/m\ < 10^^^ yi^^ we see that the 
effects ( p.24aj) on the size of the orbit is too small to be observable in satellite data. On the 
other hand the effect ( 3.24b ) leads to a mean anomaly perturbation 
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nt' + const = nt + -nt^ 2— + 2— + 3— , (3.25) 
2 \ G m m. 



corresponding to an along-track displacement^ 



2 ( G m Srris 

aoi = nat — H 1 

\ G m 2 , 



/ G m 3 m,\ n i-? i , s 

- + — + X 2.50 X lO^^cm yr'^ , 3.26 

\G m 2ms) 



which suggests that the level G/G ~ 10^^'^ yr^^ might be reachable by satellite data [as- 
suming that the ordinary sources of variation of the mass of the Earth are smaller, or, at 
least, measurable with better precision]. 

IV. CONCLUSIONS 

We have reexamined the non-Einsteinian effects in the orbital motion of low-orbit arti- 
ficial Earth satellites. We worked within the most conservative framework for alternative 
gravity theories, i.e. the assumption that gravity contains, besides the usual Einsteinian 
tensor interaction, a scalar contribution. The post-Newtonian limit of generic tensor-scalar 
gravity theories exhibit three types of non-Einsteinian effects: (i) deviations from the general 
relativistic post-Newtonian effects, parametrized by two (Eddington) parameters /3 = /3 — 1, 
and 7 = 7 — 1, (ii) combined violations of the weak and strong equivalence principles, 
parametrized by body-dependent parameters 5a = 5a + (4/3 — ^)E^^^/mAC^] and (iii) pos- 
sible slow time variabilities of Newtonian parameters, measured by G/G and rriA/TnA- 

Let us recall that the Lunar Laser Ranging (LLR) experiment has, as predicted [ll2|,p|. 



proven to be a superb tool for measuring ^Moon — ^Earth = (— 2.7±6.2) x 10~^^, i.e. for testing 
the validity of the effacement of the internal structure of the bodies in their translational 
equations of motion. The present paper indicates that the centimeter-level tracking of low- 
orbit drag-free satellites is a promising tool for measuring at the 10~^ level two independent 
combinations of (3 and 7 [/? — 27 Eq. ( |3.1| ) and (3 — ^7, Eq. ( |3.10| )] and for probing the time 
variability of G (and the masses) at the 10~^^ yr~^ level, Eq. ( |3.26| ). It is interesting to notice 
that the satellite measurement of 77 = 4/3 — 7 is complementary to the LLR measurement of 
5i — 82 in two ways: (i) conceptually it is a test of the effacement of the external universe 
in the equations of motion of a local system, and (ii) technically it gives directly access to 
the parameter r] which is "contaminated" by possible weak-equivalence-principle violations 
in the LLR observable 5^ - ^2 = - ^2 + x 4.45 x 10"^° [p^. 



9r 



The result (|3^ ) assumes that one has experimentahy access to the basic dynamical time t. 



Taking into account that time is measured by atomic clocks based on, e.g., a Bohr frequency 
cx meO? adds the terms —^7hf>/me — aja in the parenthesis on the right-hand side of Eq. ( |3.26| ). 
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A crucial difficulty in reaching the precision levels suggested by our analysis will be to 
protect the motion of the satellite from non-gravitational forcesj^. The most promising 
orbital effects, Eqs. ( |3.1| ), ( |3.10| ) and ( p.26| ), come from the build up of small effects over 
observation times of the order of one year. The drag-free system of the satellite must 
therefore be effective for correspondingly ultra low frequencies. It may well be that the 
active drag compensation systems of the type planned for the GPB satellite will fall short 
of providing the needed protection. The best system would probably be a passive drag-free 
system where the satellite floats freely within a large shell. 

If it appears impossible to use the orbital data of drag-free satellites conceived for other 
missions [like GPB or the Satellite Test of the Equivalence Principle (STEP)], it may be 
worth thinking about a dedicated mission [say the Satellite Test of the Effacement property 
(STEFF)] consisting of an optimized passive drag-free satellite: a small ball freely floating 
within a large shell which follows optically the motion of the ball and which carries some 
thrusters and GPS or DORIS receivers and/or laser corner cubes. The orbital data of such 
a mission would be of value both for geodesy and as tests of relativistic gravity. 

Other difficulties concern the knowledge of the multipole moments of the Earth which 
limit other proposals (see For instance, the present uncertainty in the Earth quadrupole 
J2, 5J2/J2 ~ 6 X 10~^ contributes an uncertainty in the predicted (Newtonian) value 
of the perigee advance which corresponds to the uninteresting level 5{(3 — 27) ~ 1. This is 
however not a serious difficulty as the high-precision satellite data that we wish to use to 
get (as by products) measurements of (3 and 7, will be primarily used anyway to refine the 
measurement of the Earth gravity field with a correspondingly high precision (see e.g. [@). 

One might think that the type of orbital tests we have been discussing will soon loose 
their interest as other planned experiments aim at reaching higher precisions in measuring 
possible deviations from General Relativity. [We have in mind, in particular. Gravity Probe 
B (which should measure 7 at the 3 x 10~^ level ||2^) and the Solar Orbit Relativity Test 
(SORT) |]2B| which aims at the 10"^ level for 7.] As these tests will not measure /?, we 



see on the contrary that orbital tests may provide our best handle on the parameter]^ (3. 
Assuming the knowledge of 7, the present analysis shows that orbital data will contain two 
independent signals, Eq. (|3.1|) and (|3.1CI|) , for measuring (5 at the lO""^ level. One can note 
that while the perigee advance of Eq. ( p.l|) is difficult to measure in the usual case of low- 
eccentricity orbits, the yearly modulation ( |3.10| ) is essentially independent of the magnitude 
of e. It remains however to see to what extent it will possible to distinguish it from the 
other yearly (gravitational) perturbations in the orbit. [Note that the phase of the signal 
( p.lOp is predicted and that the period is the anomalistic year (perihelion to perihelion).] 
Let us also remark that, contrary to many of the effects discussed in [Q, the magnitudes of 
the leading perturbations we have discussed, Eqs. ( |3.1| ), ( p.lO| ) and ( p.26| ), are independent 



One notes in particular that solar radiation pressure effects mimic exactly the yearly variation 
of Gioc but are several orders of magnitude larger. 

^^Let us note in passing that the scenario of Ref. [^ indicates that the ratio — /3/7 = /?3/4 might 
be of order 10 (if k ~ 1). 



15 



of the value of the inchnation of the orbit. A strong dependence upon / appeared only in 
the subdominant effects caused by B2, Eq. (|3.19| ). 

Finally, let us stress that the value of the approximate analytical solutions for non- 
Einsteinian perturbations that we have given above is mainly indicative of the results one 
can expect. In practice, we advise to resort to direct numerical integration of the equations of 
motion ( |2.2|) . For Agr one should use the full general relativistic post-Newtonian equations 
of motion [^, to which one can add linearly the effect of B. The accuracy with which 
we have given B in Eqs. (|2.5|) - (|2.7|) should be amply sufficient. [If one wishes one can also 
consider the non-Einsteinian effect associated with the spin angular momentum of the Earth 
S2, i.e. add the term B5, Eq. to B.] 
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